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Abstract 

We discuss the damping of tensor modes due to anisotropic stress in inflation. The effect is 
negligible in standard inflation and may be significantly large in inflation models that involve 
drastic production of free-streaming particles. 
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I. INTRODUCTION 



The spatial flatness and homogeneity of the present Universe strongly suggest that a 
period of de Sitter expansion or inflation had occurred in the early Universe pQ. During 
inflation quantum fluctuations of the inflaton field may give rise to energy density pertur- 
bations (scalar modes) |2j, which can serve as the seeds for the formation of large-scale 
structures of the Universe. In addition, a spectrum of gravitational waves (tensor modes) is 
produced from the de Sitter vacuum [3]. 

Gravitational waves are very weakly coupled to matter, so once produced they remain as 
a stochastic background till today, and thus provide a potentially important probe of the 
inflationary epoch. Detection of these primordial waves by using terrestrial wave detectors 
or the timing of millisecond pulsars [4J would indeed require an experimental sensitivity of 
several orders of magnitude beyond the current reach. However, like scalar perturbation, 
horizon-sized tensor perturbation induce large-scale temperature anisotropy of the cosmic 
microwave background (CMB) via the Sachs- Wolfe effect [5]. The recent seven-year WMAP 
anisotropy data has placed an upper limit on the contribution of tensor modes to the CMB 
anisotropy, in terms of the tensor-to-scalar ratio, that is r < 0.36 More stringent limits, 
r < 0.17 — 0.19, have been made by combining several other cosmological measurements [7]. 
In addition, the tensor modes uniquely induce CMB B-mode polarization that is the primary 
aim of on-going and future CMB experiments [8j. 

As is well-known, gravitational waves propagate freely in the expanding Universe [H]. This 
is under the assumption that the Universe is a perfect fluid. In the presence of non- vanishing 
anisotropic stress, an additional source term to the gravitational wave equation should be 
included [lO]. The effect of anisotropic stress on cosmological gravitational waves due to 
free-streaming neutrinos after the neutrino-matter decoupling in the early Universe was first 
numerically calculated in Ref. [H], and incorporated in an integro-differential equation for 
the wave propagation [12]. In fact, this equation can be also applied for any unknown free- 
streaming relativistic particles [13] . It was found that the anisotropic stress reduces the wave 
amplitude, thus lowering the tensor-mode induced CMB anisotropy and polarization [TTHT^. 

In this paper, we will discuss the effect of anisotropic stress on tensor modes in inflation. 
Here the anisotropic stress is due to free-streaming relativistic particles produced during 
inflation. The generating source of these relativistic particles could be de Sitter quantum 
fluctuations of the inflaton itself in standard slow-roll inflation [2], a thermal component 
in warm inflation [15], isolated bursts of instantaneous massless particle production [16] . 
particle production in trapped inflation in which the inflaton rolls slowly down a steep 
potential by dumping its kinetic energy into light particles at the trapping points along the 
inflaton trajectory [13118], or electromagnetic dissipation in natural inflation [T9| [20]. 

II. PARTICLE PRODUCTION IN INFLATION 

Here we assume a flat Friedmann-Robertson- Walker metric, 

ds"^ = - a^{t)dx^ = a^{ri){drf - cix^) , (1) 

where a{rj) and drj = dt/a{t) are the scale factor and conformal time respectively. For 
simplicity, we treat the inflaton energy density as approximately constant and then we 
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have 

a=-^. H' = ^P*. (2) 

We denote the energy density of the free-streaming relativistic particles x produced during 
inflation by and define a ratio, 

/ = (3) 

In the standard slow-roll inflation, a candidate for the x particle is any weakly interacting 
field whose quanta are gravitationally produced during inflation. These de Sitter quantum 
fluctuations have a characteristic energy density roughly equal to p^^ ~ H'^ [2J, so we have 
/ ~ GH^. The WMAP results [6] have set an upper limit on the inflation scale which means 
that / < 10~^°. We will see below that this small / implies too weak anisotropic stress to 
affect the tensor modes in the standard slow-roll inflation. 

However, some inflation models that involve particle production via interactions between 
and X may allow a relatively large value for /. For example, in the trapped inflation model 
with an interaction of the type (7^(0 — [E], when rolls slowly down a steep potential 
by dumping its kinetic energy into Xi particles at each trapping point at 0j along the inflaton 
trajectory, a roughly constant energy density of all Xi's, Px = J^iPxi^ maintained. This 
energy density can be estimated as p^ ~ with k^, being a characteristic energy scale given 



by fc* ~ y 9id(f)/dt), where d(j)/dt is the inflaton rolling speed. To have a successful trapped 
inflation, it is required that / < 0.1. Also, it is shown that the scattering rate of Xi particles, 
which is given by F ~ with E being the energy of Xiy is sufficiently slower than the 

expansion rate, T < H [21j. Thus, in this model the de Sitter vacuum may be populated 
with free-streaming Xi particles that generate significantly large anisotropic stress to damp 
the tensor modes. 

Another example that also provides with a constant p-^ during inflation involves an in- 
teraction, (j)/MF^^F^y, where x is a massless U{1) gauge field, F^^^ is its field strength, and 
M is a mass scale [121 120] • The growth solution for the Fourier mode of the vector potential 
with + circular polarization is found as 

1 / k 



in the interval (8^) ^ ^ k/{aH) ^ 2^, where ^ = 2{d(j) / dt) / {M H) is treated as constant. 
Hence, the energy density of the produced gauge quanta is given by 

Px = 4^ Jdke{\dA^/dr^f + e\A 
6! 



219^^2 ^3 



e2-€ . (5) 



These gauge quanta, in turn, source inflaton fluctuations which are highly nongaussian. 
The WMAP bound on nongaussianity implies that ^ ^ 3 [20]. When ^ = 3, / ~ 6.6 x 
10^ GH^ <C 1. However, the value of ^ gets to increase towards the end of inflation. If 
^ = 5.16 near the end of inflation and GH^ = 10"^*^, then we will have / ~ 0.1. The gauge 
quanta, once produced, scatter with the inflaton fluctuations with a rate given by 

r-"«S^x|^io--(-) ^, (6) 
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where the energy of the gauge particle is of order H as shown in Eq. (|4]) and the number 
density of inflaton fluctuations is approximated as < n^. As long as ^ < 5.16 and 
H < 10~^M, we reach the condition T < H, under which the gauge quanta freeze out and 
decouple from the background. 



III. GRAVITATIONAL WAVE EQUATION 

In the weak field approximation, small metric fluctuations are ripples on the background 
metric: 

gi,u = a^{ri){r]^u + h^y) , /i^^ < 1 , (7) 

where rj^j^u is the Minkowski metric and Greek indices run from to 3. In synchronous gauge, 
^00 = h^i = 0, where i runs from 1 to 3. The remaining hij contain a transverse, traceless 
tensor which corresponds to a gravitational wave or tensor mode. Henceforth we will work 
in the TT gauge, i.e., = dih"^ = and denote the two independent polarization states 
of the wave as +, x. The propagation of gravitational waves in an expanding space-time is 
well studied. In the presence of anisotropic stress, the Fourier mode equation is given by 



^ + ~^ + k%, = IQnGa'n,, , (8) 
drj^ adr] di] 

where vrjj is the Fourier mode of the TT part of the anisotropic stress tensor. 



IV. ANISOTROPIC STRESS TENSOR 

In this section, we will derive the evolution of the anisotropic stress tensor of the free- 
streaming relativistic particles x- We will follow the methodology in Ref. [10], taking into 
account the particle production in inflation and assuming that the produced x particles are 
decoupled from the background. Let A^^(x, p,t) be the phase space density of x particles, 
then the physical energy density of x particles is given by 

p^{t) = a-\t)^ J^d^^d^p\p\N^{^,p,t). (9) 

In light of the results in Sec. |Tl| in the following we will assume that the physical energy 
density is constant during inflation: dp^{t)/dt = 0. This suggests that we should deal with 
a re-scaled phase space density instead, defined by 

n^{^,p,t)=a-\t)N^{^,p,t). (10) 

In the absence of collisions, the re-scaled phase space density satisfies a Boltzmann equa- 
tion in a metric 5'ij(x, t), 

dn^ ^ dn^jf_ ^ dn^ p>p^ dgjk _ ^ 
dt (9x* dpi 

where = g'^^Pj and p° = \f(f^^PiPj- We now consider small perturbation, 

gij = a^5ij + 5gij, = n^{ap^) + 5n^, (12) 
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(11) 



where with p = |p| is just the re-scaled local phase space density, which has a well- 

defined energy spectrum as exemplified in Eq. ([5]). To first order in metric and density 
perturbation, Eq. (11) reads 



dn^{p) 



dt 
d6n^ 
~dt 



+ 



= 0, 

Pi d5n^ 
a dx^ 



p dn^jp) ^ . d / 6gi 
2 dp ^'^'dt [ a2 



(13) 
(14) 



where pi = Pi/p- 

Using Eq. ([T]) for the metric perturbation, we write down the tensor component of Eq. (14) 

as 



d5n 



pdn^{p) ^ ^ dh, 



dt ' a dx^ 2 dp ' dt ' 
Let us introduce a dimensionless re-scaled intensity perturbation defined by 



Then, Eq. (15) becomes 



p^J{x,p,t) = J dpA7ip^5n^{x,p,t), 
Px = J dpAnp^n^ip). 



(15) 

(16) 
(17) 

(18) 



where we have used integration by part and assumed that n^iS^) = n^{oo) = 0. We can then 
construct the spatial component of the stress tensor perturbation as 



6T 



XJ 



d^p6n^{:si,p,t)ppiPj, 



(19) 



which contributes to the anisotropic stress tensor tt^. Following the same steps in Ref. [TU] . 
the free-streaming solution for the anisotropic stress tensor of x particles in the presence of 
gravitational waves is found as 



TT, 



where rji is some initial time and the kernel is given by 



K(u) 



smii 



3cosM Ssin-u 

^ + ^ 



(20) 



(21) 



The integro-differential equation ([8j) with Tiij given by Eq. (20) has been solved for the case 
in which x particles are free-streaming neutrinos in the early Universe [T2I4H] . One would 



anticipate that the free-streaming solution of the anisotropic stress tensor (20) is a back 
reaction to the wave equation and thus reducing the wave amplitude. 

In trapped infiation [T7j or axionic infiation with 10^ U{1) gauge fields [11], the time-scale 
of particle production is much shorter than the expansion time, H~^. For instance, trapped 
infiation produces particles in a time-scale of order <^ H~^. Therefore, infiation begins 
very shortly after particles are copiously produced. Let rji be the moment when infiation 
begins. The initial condition, ^^^{rii) = 0, that we have assumed in Eq. (20) is then justified. 
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Note that generically anisotropic stress should exist before inflation, i.e. Tiijirji) ^ 0, due 
to the fact that we do not really know the initial condition for inflation since here we do not 
have a physical model before inflation. However, soon after inflation starts, this pre-existing 
anisotropic stress has decayed and become vanishingly small. Since we are mainly interested 
in X particles produced during inflation, we have assumed that the generic anisotropic stress 
is absent at the beginning of inflation, namely 7Tij{r]i) = 0. Otherwise, we will need to 
consider the effect of this generic anisotropic stress on gravitational waves in a brief period 
after the start of inflation. 



V. DAMPING IN INFLATION 



Let us decompose 

hijiv,^) = {87iGy^Hk-h{7],k)e,,{k;X), (22) 

where ejj(k; A) is the polarization tensor with A = +, x. Here we have introduced a dimen- 
sionless wave amplitude h that is assumed to be the same for both polarizations. Then, 
Eq. ^ becomes 

d'h 2dh . 24/ r ,.dh 

K{u — u)—du. (23) 

(Jj LL 



du^ u du 



+ h 



where / = Pxl is a constant. The homogeneous solution of Eq. (23) is known. Selecting 
the Bunch-Davis vacuum [22], it is given by 



/lo = —(1 + "I'U) e 



(24) 



At the end of inflation (i.e., u — )■ 0), |/io| — ?■ 1- This reproduces the scale-invariant power 
spectrum predicted in standard slow-roll inflation. 

Since the damping effect is expected to be secondary, we can make the Born approxi- 



mation to replace dh/du' in the damping term in Eq. (23) by dh^/du' . Then, we use the 

(25) 



retarded Green's function method to find the particular solution 

K = 



^ dii' /■"' dhn 

24// ^G{u-u') du"K{u'-u") ' 



du" 

where the Green's function is constructed from the homogeneous solution, given by 

1 \ , u — w 



G{u - u') 



u 
u' 



uw 



smm — u 



cos Ui — u 



uw 



(26) 



Now we set a = 1 at the time when inflation starts. This fixes ui = —k/H with k/H = 1 
corresponding to the length scale that crosses the horizon at the onset of inflation. We then 



numerically evaluate the integral (25) by letting m — )■ to obtain hp for a range of values of 
k/H. In Fig. [l| assuming that 24/ = 1, we plot the total wave amplitude, \h\ = \hQ + hp\, 
at the end of inflation {u — )■ 0) against the wave number, k/H. The figure shows that the 
wave amphtude is reduced and becomes asymptotically flat at large k. 

Apparently, the wave amplitude is enhanced for k/H < 3. However, we do not expect 
that the present work would give precise results for these low fc-modes. The enhancement 
should be an artifact due to the use of the approximation, dh/du' ~ dho/du', and the choice 
of the Bunch-Davis vacuum near the initial time, Wj. A quick way to resolve this problem. 
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albeit artificial, is to integrate the damping effect on each fc-mode for the time interval from 
the start of inflation to the horizon crossing time, which means by setting m = — 1 in Eq. (25 ). 
As such, = 1 for k/H = 1 by default. Some improvements can also be made such as 
using an advanced solving method for the integro-differential equation ^ and considering 
a smooth transit to the inflationary phase. 




FIG. 1: The wave amplitude \h\ versus the wave number k/H at the end of inflation by the use of 
the Born approximation in Eq. (25), where 24/ = 1 is assumed and k/H = 1 corresponds to the 



length scale that leaves the horizon at the onset of inflation. 



VI. CONCLUSION 



We have discussed the damping effect of anisotropic stress on tensor modes due to free- 
streaming relativistic particles produced during inflation. The damping increases with the 
ratio of the particle energy density and the de Sitter vacuum energy, which ranges from 
10~^° in standard inflation to about 0.1 in inflation models that involve drastic particle 
production such as the trapped inflation or axionic inflation. In these inflation models, the 
particle production may significantly reduce the amplitude of the tensor modes. 

Recently new sources of gravitational waves during inflation have been proposed. They 
are anisotropic stress induced by quantum energy stress of conformal fields and their 
associated fluctuations |21] , by the Bremsstrahlung from the particle production events , 
and by the produced gauge field quanta that couple to inflaton All of these produced 

gravitational waves may experience the damping effect considered in the present work if 
copious free-streaming relativistic particles are also produced during inflation. As such, to 
properly take into account the damping, one needs to consider the full integro-differential 
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equation, 



(Phii 2dadhij ,9- ^ 9 P 



'■J 



(27) 



where the damping term is taken from Eq. (20) and nfj''" denotes a new source term for 
the anisotropic stress. When / <^ 1, the dampmg term can be neglected and the equation 



reduces to that considered in Refs. [2n [231 - I26] . Otherwise, Eq. (27) should be solved self- 
consistently to obtain the damped tensor power spectrum. 



At last we note that in Ref. [26] 
they have considered the production of gravitational waves at the interferometric scales 
during the final stage of inflation when ^ ~ 5 — 6. This large value of ^ may imply that 
/ ~ 1, thus resulting in a large damping on the tensor modes. 
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